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Abstract 


This  report  addresses  the  problem  of  discrimination  power  in  target  tracking  applications. 
More  specifically,  a  closed-loop  approach  to  adapt  the  sensing  and  tracking  operations  is 
proposed  and  compared  to  the  conventional  open-loop  and  static  approach.  The  objective 
is  to  control  and  maintain,  over  a  certain  volume  of  interest  and  by  way  of  clustering  and 
scheduling  strategies,  the  level  of  discrimination  power  required  by  the  mission  objectives. 
The  control  strategy  is  based  on  two  cascade  loops.  The  outer  loop  uses  clustering  tech¬ 
niques  to  characterize  the  volume  of  interest  in  terms  of  discrimination  power.  This  high 
level  information  is  exploited  by  the  inner  loop  to  compute  optimal  track  update  and  sen¬ 
sor  scheduling  strategy.  The  presented  results  show  that  the  discrimination  power  can  be 
improved  by  adapting  the  target  tracking  operations.  This  improvement  could  benefit  tac¬ 
tical  military  surveillance  operations  such  as  contact/track  correlation,  target  engagement, 
target  identification  and  classification. 


Resume 


Le  travail  presente  dans  ce  rapport  porte  sur  le  probleme  de  la  capacite  discriminatoire  dans 
des  applications  de  pistage  de  cibles.  Plus  precisement,  une  approche  en  boucle  ferrnee,  qui 
perrnet  d’adapter  les  operations  de  detection  et  de  pistage,  est  proposee  et  comparee  a  une 
methode  conventionnelle  statique  en  boucle  ouverte.  L’objectif  poursuivi  est  de  controler  et 
de  maintenir,  pour  un  certain  espace  d’interet  et  a  l’aide  de  techniques  de  regroupement  de 
donnees  et  de  gestion  de  taches,  le  niveau  de  capacite  discriminatoire  requis  pour  repondre 
aux  objectifs  d’une  mission.  La  strategie  de  controle  est  basee  sur  deux  boucles  imbriquees. 
La  boucle  externe  utilise  des  techniques  de  regroupement  de  donnees  pour  caracteriser 
l’espace  d’interet  sur  la  base  de  la  capacite  discriminatoire.  Cette  information  de  haut 
niveau  est  alors  exploitee  dans  une  boucle  interne  qui  determine  une  strategie  optimale 
de  gestion  de  la  rnise  a  jour  des  pistes.  Les  resultats  obtenus  montrent  que  la  capacite 
discriminatoire  peut  etre  amelioree  en  adaptant  le  pistage  de  cibles.  L’ amelioration  ainsi 
apportee  pourrait  profite  r  aux  operations  militaires  de  surveillance  telles  que  la  correlation 
mesures/ pistes,  l’identification,  la  classification  et  l’engagement  de  cibles. 
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Executive  summary 


Discrimination  power  control  for  adaptive  target  tracking 
applications 

A.  Benaskeur,  F.  Rheaume;  DRDC  ValcartierTR  2008-016;  Defence  R&D  Canada  - 

Valcartier;  July  2008. 

In  military  Command  and  Control  applications,  surveillance  and  target  tracking  aim  at  pro¬ 
viding  accurate  and  timely  identification,  classification,  and  kinematics  information  about 
the  entities  within  the  volume  of  interest.  Modern  surveillance  and  tracking  systems  have 
been  making  increasing  use  of  data  fusion  technology  and  tools  to  achieve  their  mission’s 
goals  and  objectives.  This  work  addresses  the  problem  of  discrimination  power  within  the 
framework  of  adaptive  data  fusion. 

In  this  work  discrimination  power  is  defined  as  the  capability  to  distinguish  objects  based 
on  their  track  information.  Two  different  levels  of  discrimination  power  are  defined,  the 
signal  and  the  statistical  level.  The  presented  target  tracking  problem  concentrates  on 
the  statistical  level.  According  to  that,  different  discrimination  power  measures  are  sug¬ 
gested:  coverage  intervals,  Mahalanobis  distance  and  chi-squared  testing.  Because  the 
discrimination  power  is  understood  as  a  critical  factor  for  operations  such  as  contact /track 
correlation,  target  identification  and  classification,  we  propose  a  method  that  adaptively 
controls  it  within  the  tracking  system.  This  discrimination  power  control  should  also  bene¬ 
fit  to  target  engagement  operations,  which  require  a  high  discrimination  power  in  order  to 
maximize  the  chance  of  threat  neutralization  and  minimize  the  risk  of  collateral  damages. 

The  suggested  discrimination  power  control  method  rely  on  concepts  of  adaptive  data  fu¬ 
sion.  The  adaptation  uses  a  two-level  cascade  loop,  the  objective  of  which  is  controlling 
discrimination  power  over  a  certain  volume  of  interest.  The  outer-loop  characterizes  the 
whole  volume  of  interest  based  on  the  concept  of  discrimination  power.  Using  clustering 
and  classification  techniques,  distinct  regions  are  created  and  managed  dynamically  based 
on  the  targets’  spatial  distribution.  The  inner-loop  controls  discrimination  power  over  the 
volume  of  interest  by  an  appropriate  selection  of  the  track  update  frequency  and  sensor 
scheduling  within  each  distinct  region.  Three  different  scheduling  methods  are  proposed: 
time-slice,  minimum  intra-cluster  discrimination  power  and  round-robin.  All  the  methods 
rely  on  clusters  of  tracks,  but  only  time-slice  and  minimum  intra-cluster  discrimination 
power  adaptively  control  the  track  update  periods.  Round-robin  uses  a  classical  open-loop 
strategy. 

A  set  of  scenarios,  with  varying  target  densities,  was  simulated  to  test  different  discrimi¬ 
nation  power  control  strategies.  The  scenarios  consider  a  single  sensor  that  can  switch  the 
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direction  of  its  beam  instantly.  The  presented  results  have  shown  that  the  discrimination 
power  can  be  improved  by  adapting  the  target  tracking  operations.  The  time-slice  schedul¬ 
ing  algorithm  proved  to  be  better  than  minimum  intra-cluster  discrimination  power  and 
than  the  open-loop  round-robin  strategy.  The  time-slice  approach  was  able  to  maintain  a 
higher  discrimination  power  in  regions  with  high  target  densities. 
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Dans  les  applications  militaires  de  commandement  et  controle,  les  operations  de  surveillance 
et  de  pistage  de  cibles  visent  a  fournir  de  maniere  precise  et  rapide  de  l’information  sur 
l’identification,  la  classification  ainsi  que  sur  la  cinematique  des  cibles  presentes  dans  un 
espace  d’interet.  Les  systemes  de  surveillance  modernes  font  de  plus  en  plus  appel  aux 
technologies  et  aux  outils  de  fusion  de  donnees  pour  realiser  les  objectifs  propres  a  leur 
mission.  Ce  travail  porte  sur  le  probleme  de  la  capacite  de  discernement,  presente  dans  le 
cadre  de  la  fusion  adaptative  de  donnees. 

La  capacite  de  discernement  est  definie,  dans  ce  travail,  comrne  la  capacite  de  distinguer  des 
objets  entre  eux,  sur  la  base  de  l’infornraiton  contenue  dans  les  pistes.  La  capacite  de  dis¬ 
cernement  se  divise  en  deux  niveaux,  le  niveau  du  signal  et  le  niveau  statistique.  Ce  dernier 
sera  etudie  dans  ce  travail.  A  ce  propos,  differentes  mesures  de  la  capacite  de  discernement 
sont  suggerees  :  intervalles  de  couverture,  distance  de  Mahalanobis  et  test  du  chi-carre.  La 
capacite  de  discernement  etant  consideree  comrne  un  facteur  critique  pour  les  operations 
telles  que  la  correlation  mesure/piste,  l’identification  et  la  classification,  une  methode  pour 
controler  adaptativement  la  capacite  de  discernement,  dans  un  systeme  de  pistage,  est  pro- 
posee.  Le  controle  adaptatif  de  la  capacite  de  discernement  pourrait  egalement  profiter  aux 
operations  d’engagement  de  cibles.  Ces  dernieres  requierent  un  grande  capacite  de  discer¬ 
nement  afin  de  maximiser  la  probability  de  neutralisation  de  la  menace  et  de  minimiser  le 
risque  de  dommages  collateraux. 

La  methode  de  controle  de  la  capacite  de  discernement  proposee  se  base  sur  des  concepts 
de  fusion  adaptative  de  donnees.  L’adaptation  est  realisee  grace  a  deux  boucles  en  cascade. 
L’objectif  poursuivi  est  de  controler  la  capacite  de  discernement  pour  un  certain  volume 
d’interet.  La  boucle  externe  definit  le  volume  d’interet  sur  la  base  de  la  capacite  discrimi- 
natoire.  A  l’aide  de  techniques  de  regroupement,  differentes  regions  de  l’espace  sont  creees 
et  gerees  dynamiquement  selon  la  distribution  spatiale  des  cibles  observees.  La  boucle  in¬ 
terne  controle  la  capacite  discriminatoire  par  la  determination  des  strategies  de  rnise  a  jour 
des  pistes  pour  chaque  region  identifies  du  volume  d’interet.  Trois  differentes  strategies  de 
rnise  a  jour  sont  proposees  :  “tranche  de  temps” ,  “capacite  de  discernement  minimale  intra- 
groupe”  et  “a  tour  de  role”.  Toutes  ces  methodes  se  basent  sur  le  regroupement  de  cibles 
rnais  seulement  “tranche  de  temps”  et  “capacite  de  discernement  minimale  intra-groupe” 
sont  adaptatives.  La  methode  “a  tour  de  role”  utilise  la  strategie  classique  en  boucle  ouverte. 
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Des  scenarios  avec  differentes  densites  de  cibles  ont  ete  simules  afin  de  tester  les  strategies 
developpees  sur  le  controle  de  la  capacite  discriminatoire.  Un  capteur  unique,  capable  de 
diriger  son  faisceau  de  nraniere  instantanee,  est  utilise.  Les  resultats  obtenus  nrontrent  qu’il 
est  possible  d’anreliorer  la  capacite  discriminatoire  en  adaptant  les  operations  de  pistage. 
L’algorithnre  “tranche  de  temps”  s’est  revele  meilleur  que  celui  base  la  “capacite  de  discer- 
nenrent  mininrale  intra-groupe”  et  que  l’algorithme  “a  tour  de  role”.  L’algorithnre  “tranche 
de  temps”  a  ete  en  rnesure  de  nraintenir  une  capacite  de  discernenrent  plus  elevee  dans  les 
regions  a  haute  densite  de  cibles. 


VI 


DRDC  Valcartier  TR  2008-016 


Table  of  contents 


Abstract .  i 

Resume .  i 

Executive  summary  .  iii 

Sommaire .  v 

Table  of  contents .  vii 

List  of  figures .  ix 

List  of  tables .  xi 

1  Introduction .  1 

2  Data  fusion  in  military  operations  .  3 

2.1  Adaptive  data  fusion .  5 

2.2  Target  Tracking .  6 

2.2.1  Kalman  filtering .  7 

2.2.2  Adaptation  in  tracking .  8 

3  Discrimination  power .  9 

3.1  Signal  level  discrimination  power  .  10 

3.2  Statistical  level  discrimination  power .  11 

3.2.1  Discrimination  power  metrics .  12 

3. 2. 1.1  Probability .  12 

3. 2. 1.2  Coverage  intervals .  13 

3. 2. 1.3  Mahalanobis  distance .  15 

3. 2. 1.4  x2  test  for  assessing  discrimination  power  .  16 

DRDC  Valcartier  TR  2008-01 6  vii 


4  Discrimination  power  control 


17 


4.1  Simple  track  update  algorithm .  17 

4.2  Clustering-based  control  method .  17 

4.2.1  Hierarchical  clustering  of  tracks .  18 

4.2.2  Sensor  scheduling  heuristics .  20 

4. 2. 2.1  Time-slice .  21 

4. 2. 2. 2  Minimum  Intra-Cluster  Discrimination  Power  (MICDP)  22 

4. 2. 2. 3  Round-robin .  22 

5  Results  and  Discussion .  23 

5.1  Sensor  Model  .  23 

5.2  Performance  measures  .  23 

5.2.1  Discrimination  Conflict  .  24 

5.2.2  Average  track-to-track  Mahalanobis  distance .  24 

5.2.3  Average  error  covariance  matrix  determinant .  24 

5.3  Scenarios .  24 

5.3.1  Scenario  A  .  26 

5.3.2  Scenario  B  .  26 

5.3.3  Scenario  C  .  28 

5.4  Notes  on  clustering  and  scheduling  .  31 

6  Conclusion  .  35 

List  of  symbols .  37 

List  of  acronyms .  39 

Glossary .  41 

References .  43 


viii 


DRDC  Valcartier  TR  2008-016 


List  of  figures 


Figure  1:  JDL  model  of  data  fusion .  3 

Figure  2:  Open-loop  data  collection  and  fusion  (no  adaptation) .  4 

Figure  3:  Closed-loop  data  collection  and  fusion  (with  adaptation) .  5 

Figure  4:  Two-level  Cascade  Control  Loop  .  7 

Figure  5:  Sensor  with  accuracy  a  (Gaussian  probability  distribution)  and 

resolution  <p .  9 

Figure  6:  Discrimination  power  levels .  10 

Figure  7:  Two  targets  observed  in  a  ID  space .  11 

Figure  8:  Discrimination  power  for  two  state  estimates  .  11 

Figure  9:  95%  coverage  interval  (Cl)  of  an  estimate  x .  13 

Figure  10:  Confidence  intervals  and  discrimination  power  -  CI\  and  CI2  do  not 

overlap,  the  discrimination  power  is  fair .  14 

Figure  11:  Confidence  intervals  and  discrimination  power  -  CI\  and  CI2  overlap, 

the  discrimination  power  is  too  low .  14 

Figure  12:  Statistical  distance  d  as  a  track  dissimilarity  metric .  15 

Figure  13:  Clustering/discrimination  power-based  Two-Level  Control  Loop .  18 

Figure  14:  Distance  tree  and  hierarchical  formation  of  clusters .  19 

Figure  15:  Track  Quality  as  function  of  scheduling  strategy  for  the  same  number  of 

updates .  22 

Figure  16:  Trajectories  of  the  4  targets  in  Scenario  A  (‘X’  is  the  final  position)  ...  26 

Figure  17:  Mahalanobis  distance  between  tracks  of  targets  1  and  2.  The  ground 

truth  distance  for  targets  1  and  2  is  also  illustrated .  27 

Figure  18:  Track  quality  yielded  by  round-robin  scheduling  strategy .  28 


DRDC  Valcartier  TR  2008-016 


IX 


Figure  19:  Track  quality  yielded  by  the  time-slice  discrimination  power 

control-based  scheduling  strategy .  29 

Figure  20:  Sample  of  trajectories  of  the  fifteen  targets  in  Scenario  C  (‘X’  is  the 

final  position) .  30 

Figure  21:  Number  of  clusters  in  terms  of  time .  31 

Figure  22:  Average  track-to-track  Mahalanobis  distance  (AM) .  32 

Figure  23:  Total  number  of  Discrimination  Conflicts  (DCu  ..,15]) .  33 

Figure  24:  Average  determinant  of  error  covariance  matrix  ( ADET )  33 


x 


DRDC  Valcartier  TR  2008-016 


List  of  tables 


Table  1:  Closest  pairs  of  targets  in  scenario  B .  27 

DRDC  Valcartier  TR  2008-01 6  xi 


This  page  intentionally  left  blank. 


XII 


DRDC  Valcartier  TR  2008-016 


1  Introduction 


In  military  Command  and  Control  (C2)  applications,  surveillance  and  target  tracking  [1, 
2,  3,  4]  aim  at  providing  accurate  and  timely  detection,  identification,  classification,  and 
kinematics  information  about  the  entities  within  the  Volume  Of  Interest  (VOI).  Modern 
surveillance  and  tracking  systems  have  been  making  increasing  use  of  data  fusion  tech¬ 
nology  and  tools  [5,  6,  7]  to  achieve  their  mission’s  goals  and  objectives.  To  gain  higher 
performance,  a  modern  data  fusion  system  needs  an  active  feedback  or  adaptation.  The 
adaptation  may  concern  the  data  fusion  process  itself  or  the  related  sensor  management 
problem.  This  corresponds  to  Level  4  of  the  JDL  model  proposed  by  the  US  DoD  Joint 
Directors  of  Laboratories  sub-panel  [8] .  Level  4  has  been  removed  from  the  latest  version  [9] 
of  the  JDL  model,  since  it  is  viewed  as  a  resource  management  problem. 

This  work  addresses  the  problem  of  discrimination  power  characterization  and  control 
within  the  broader  paradigm  of  adaptive  data  fusion  [10,  11].  More  specifically,  the  report 
presents  the  problem  of  adapting  the  target  tracking  operations  and  attempts  to  demon¬ 
strate  the  benefits  of  such  adaptation  compared  to  the  classical  open-loop  operation  mode. 
Note  that  our  definition  of  adaptive  data  fusion  corresponds  of  the  definition  of  Level  4  in 
the  two  previous  versions  of  the  JDL  [8,  12]  and  to  the  concept  of  levels’  refinement  in  the 
current  version  [9]. 

Target  tracking  provides  estimations  of  the  targets’  states  represented  by  probability  dis¬ 
tributions  of  the  states,  which  reveal  the  uncertainty  about  the  state  values.  In  situations 
involving  multiple  targets,  the  latter  may  come  too  close  to  be  clearly  distinguishable  from 
each  other  by  the  tracking  system.  Targets  are  said  to  be  distinguishable  when  the  overlap 
between  their  spatial  probability  distributions  is  below  a  specific  level.  The  non-overlapping 
portion  of  the  spatial  probability  distributions  expresses  a  dissimilarity  relation  measured  in 
terms  of  the  distance  between  the  spatial  probability  distributions  [13].  This  dissimilarity 
relation  among  tracks  is  referred  to,  in  this  report,  as  discrimination  power.  The  latter 
is  understood  as  a  critical  factor  for  operations  such  as  contact/track  correlation,  target 
identification  and  classification.  Its  impact  may  even  be  more  critical  for  target  engagement 
operations.  Given  this  criticality,  a  high  discrimination  power  is  often  required  to  maximize 
the  chance  of  threat  neutralization  and  minimize  the  risk  of  collateral  damages. 

The  above  mentioned  discrimination  power  notion  is  used  as  an  adaptation  enabler  for  tar¬ 
get  tracking  operations.  This  adaptation  is  achieved  thanks  to  a  two-level  cascade  loop,  the 
objective  of  which  is  controlling  discrimination  power  over  a  certain  volume  of  interest.  The 
outer-loop  characterizes  the  whole  volume  of  interest  based  on  the  concept  of  discrimina¬ 
tion  power.  Using  clustering  and  classification  techniques,  distinct  regions  are  created  and 
managed  dynamically  based  on  the  targets’  spatial  distribution.  The  inner-loop  controls 
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discrimination  power  over  the  volume  of  interest  by  an  appropriate  selection  of  the  track 
update  frequency  and  sensor  scheduling  within  each  distinct  region. 

The  remaining  part  of  the  report  is  organized  as  follows.  First,  Chapter  2  introduces  data 
fusion  in  military  operations.  Chapter  3  presents  and  defines  the  discrimination  power 
problem.  Thereafter,  Chapter  4  discusses  the  control  of  discrimination  power  in  target 
tracking  as  a  special  application  of  the  adaptive  data  fusion  paradigm.  Finally,  target 
tracking  scenarios  are  presented  in  Chapter  5,  along  with  the  results  of  the  discrimination 
power  control  strategies  that  are  developed  in  Chapter  4. 


2 
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2  Data  fusion  in  military  operations 


By  reducing  uncertainty  in  the  existing  pieces  of  information  and  providing  means  to  infer 
about  the  missing  pieces,  data  fusion  supports  the  operators  in  compiling  and  analyzing  the 
tactical  picture  and  ultimately  improving  the  situation  awareness  of  the  decision  makers. 


According  to  Llinas  et  al.  [9],  in  their  revision  of  the  model  proposed  initially  by  the  US 
DoD  Joint  Directors  of  Laboratories  (JDL)  sub-panel  [8],  the  data  fusion  process  is  subdi¬ 
vided  into  five  levels  where  each  succeeding  level  deals  with  a  higher  level  of  information 
abstraction  (see  Figure  1). 


Level  0  -  is  the  sub-object  assessment  level  and  deals  with  the  estimation  and  prediction 
of  signal/object  observable  states  on  the  basis  of  pixel/signal  level  data  association 
and  characterization.  This  includes  the  signal  detection  and  feature  extraction. 

Level  1  -  also  referred  to  as  the  object  assessment  level,  uses  the  sensor  data  from  Level  0 
to  optimally  estimate  the  current  kinematic  properties  of  the  target,  and  predict  their 
future.  It  also  makes  inferences  about  the  target’s  identity  and  other  key  attributes. 
The  output  of  Level  1  is  an  aggregated  composite  tactical  picture. 

Level  2  -  concerns  the  situation  assessment  issue  and  leads  to  a  more  symbolic  represen¬ 
tation  of  the  environment  and  the  relationships  among  the  key  entities  and  the  events 
in  it. 

Level  3  -  At  the  highest  level  of  data  fusion  is  the  impact  assessment  level,  or  Level  3, 
that  projects  the  current  situation  into  the  future  and  infers  about  the  impact  of  the 
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assessed  situation,  the  vulnerability  and  the  force  capabilities. 


The  data  fusion  process  had  often  been  portrayed  as  a  purely  passive  and  open  loop  treat¬ 
ment  that  simply  transforms  the  pieces  of  information  it  receives  (see  Figure  2).  The  infor¬ 
mation  gathered  by  the  sensing  resources  undergoes  a  fusion  process  without  any  feedback 
to  the  sensors  (sensor  management)  or  to  the  processing  (fusion  adaptation). 

Figure  2  illustrates  the  fusion  process  without  feedback  (i.e.,  open  loop).  Normal  operating 
procedures  control  the  sensors  and  the  data  they  collect  is  used  in  the  fusion  process  to 
arrive  at  a  high  level  analysis  of  the  situation  of  interest.  There  is  no  formal  provision  for 
adjusting  the  sensing/processing  process  to  achieve  higher  performance. 


Figure  2:  Open-loop  data  collection  and  fusion  (no  adaptation) 


However,  the  effectiveness  of  an  advanced  system  is  not  only  determined  by  the  capabilities 
of  its  individual  functionalities  and  resources  alone,  but  also  by  the  effectiveness  of  the 
whole  system  integration.  This  integration  must  focus  on  a  cooperative,  synergistic,  and 
efficient  utilization  of  all  of  the  available  resources/processing  components.  Therefore,  to 
gain  higher  performance,  a  modern  surveillance  system  needs  many  additional  functions, 
and  the  most  essential  is  an  active  feedback  or  adaptation  (see  Figure  3). 
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Figure  3:  Closed-loop  data  collection  and  fusion  (with  adaptation) 


Figure  3  illustrates  the  fusion  process  when  adaptation  and  sensor  management  are  used  in 
a  feedback-based  sensing  and  processing  strategy.  To  make  adjustments  to  the  sensing  and 
processing  operations,  adaptation  and  sensor  management  use  outputs  from  higher  fusion 
processes.  This  corresponds  to  the  definition  and  the  role  of  Level  4  of  the  JDL  model.  It 
has  to  do  with  how  to  best  manage,  coordinate  and  organize  the  processing  at  lower  levels 
and  the  sensing  resources  in  a  manner  that  synergistically  improves  the  data  fusion  [14,  15]. 
Level  4  defines  the  adaptive  data  fusion  and  the  related  sensor  management  problems  that 
are  treated  in  the  sequel. 

2.1  Adaptive  data  fusion 

Generally  speaking,  a  system  is  called  adaptive  if  it  can  analyze  its  own  performance  and 
dynamically  reconfigure  itself  to  compensate  for  changes  in  its  context.  The  changes  might 
be  from  the  environment  it  evolves  in,  its  objectives  and/or  its  requirements.  Instead  of 
being  developed  for  a  specific  situation,  adaptive  systems  are  therefore  able  to  handle  a  wide 
class  of  situations  defined  by  a  set  of  structural  constraints  on  their  context.  Adaptation 
in  data  fusion  has  been  the  object  of  a  growing  interest  during  the  last  few  years  and  is 
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still  the  least  mature  of  the  JDL  levels  [16,  17,  10].  It  is  expected  that  within  the  data 
fusion  context,  the  introduction  of  feedback  and  adaptation  concepts  would  result  in  an 
increased  performance.  Still,  this  functionality,  which  would  make  all  the  difference,  is 
often  missing  in  modern  data  fusion  systems.  It  is  worth  noting  that  sensor  management 
represents  a  major  component  of  the  adaptation  data  fusion  problem  on  which  most  of  the 
recent  research  effort  has  been  focused  [18,  14,  19]. 

An  adaptive  system  needs  to  perceive  the  environment  it  operates  in  and  uses  this  knowledge 
to  produce  appropriate  actions  to  achieve  its  goals  [20] .  The  system  must  take  into  account 
the  possible  dynamic  nature  of  the  environment  it  interacts  with  in  order  to  reach  its  goals 
throughout  a  wide  range  of  situations,  with  the  desired  level  of  performance. 

Therefore,  handling  the  adaptation  problem,  in  data  fusion  applications,  presumes  that  the 
required  performance  of  the  data  fusion  system  can  be  specified  quantitatively  to  define  the 
adaptation  goals  and  objectives.  A  metric  can  then  be  calculated  and  used  to  evaluate  the 
data  fusion  performance  ( i.e .,  the  deviation  from  the  desired  behavior).  This  allows  the 
adaptive  data  fusion  system  to  monitor  its  behavior  and  detect  deviations  from  commit¬ 
ments  or  new  opportunities.  Actions  can  then  be  taken  to  make  the  system  attain  its  goals 
with  the  desired  level  of  performance. 

As  illustrated  by  Figure  4,  the  performance  measurement  and  control  process  can  be  per¬ 
formed  recursively  at  different  levels  of  abstraction,  where  the  loop  of  level  n  +  1  sets  the 
objectives  for  the  loop  of  level  n.  The  latter  selects  the  appropriate  actions  to  achieve  those 
objectives.  Associated  with  each  adaptation  loop  is  a  performance  measure  that  provides 
the  necessary  feedback  from  the  environment. 

2.2  Target  Tracking 

Tracking  operations,  as  part  of  the  situation  analysis  process,  is  aimed  at  providing  accu¬ 
rate  and  timely  kinematics,  identification,  and  classification  information  about  the  entities 
within  the  VOI.  This  information  is  referred  to  as  the  compiled  Tactical  Picture  (TP).  The 
military  typically  operates  in  demanding,  dynamic,  semi-structured  and  large-scale  envi¬ 
ronments.  This  reality  makes  it  difficult  to  detect,  recognize/classify  and  track  accurately 
all  entities  within  the  VOI,  thus  diminishing  the  capabilities  to  react  and  response  properly 
to  the  ones  that  pose  actual  threat.  This  can  be  very  critical  to  own- force  survival  in  high 
target-density  operations,  such  as  in  littoral  environments. 

In  this  report,  target  tracking  concentrates  on  kinematics1,  where  the  sensor  data  is  used  to 
optimally  estimate  the  current  kinematical  properties  of  the  target  and  predict  their  future 

1  Target  identification  and  classification  is  not  part  of  this  work. 
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Figure  4:  Two-level  Cascade  Control  Loop 


positions.  A  widely  used  algorithm  for  estimating  time-varying  target  states  is  the  Kalman 
filter,  which  will  be  described  below. 

2.2.1  Kalman  filtering 

For  time-varying  quantities  and  for  linear  and  Gaussian-noise  systems,  the  Kalman  filter 
provides  optimal  state  estimation  [21,  3].  However,  in  most  applications,  there  may  be 
some  non-linearity  that  can  be  associated  with  the  process  model,  the  observation  model  or 
with  both.  In  such  cases,  a  modified  version  of  the  Kalman  Filter  known  as  the  Extended 
Kalman  Filter  (EKF)  must  be  used  [22].  With  the  EKF,  the  non-linear  functions  are 
linearized  approximately  using  partial  derivatives  of  the  non-linear  functions,  so  that  the 
state  transition  and  observation  models  need  not  necessarily  be  linear  functions  of  the 
state.  The  following  give  the  equations  of  the  dynamical  model  and  of  the  EKF  used  by  the 
underlying  tracking  algorithm.  The  discrete-time  dynamical  model  of  the  targets  is  given 
by 


xk+1  =  f(xk)  +  Tvk  (1) 

Zk+ 1  —  h(xk+ 1)  T  wk+i  (2) 

where  x  is  the  n-dimension  state  vector,  v  is  the  process  noise  with  covariance  matrix  Q,  z 
is  the  measurement  vector  and  w  is  the  measurement  noise,  whose  covariance  matrix  will  be 
denoted  R.  T  is  the  discretized  continuous  time  process  noise  transition  matrix.  Thus  the 
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EKF  will  account  for  both  the  uncertainties  resulting  from  the  process  and  measurement 
noises.  The  time-update  equations  for  the  EKF  algorithm  are  given  by 

•^k+l\k  Fk/^'k  (3) 

Pk+i\k  =  FkPkFk  +  TQr7  (4) 

where  F \  is  the  Jacobian  of  /.  The  above-given  predicted  values  are  used  to  calculate  the 
updated  version  of  the  state  and  the  corresponding  error  covariance  matrix. 

-^fe+i|fc+i  =  P k+l\k  +  (5) 

^>k+l\k+l^k+1\k+1  =  ^k+l\k^k+l\k  +  zk+l  (6) 

with 

P-1  =  Hl+lR-lHk+x  (7) 

and  where  Hk  is  the  Jacobian  of  hk ■  k  represents  the  discretized  time  instant. 

The  next  chapters  will  show  how  the  tracking  process  can  be  adapted  to  control  the  dis¬ 
crimination  power  of  the  tracked  targets. 

2.2.2  Adaptation  in  tracking 

Adaptation  in  the  specific  context  of  target  tracking  aims  at  producing  a  system  that  can 
readily  adapt  to  changing  operating  environment  and  needs.  An  adaptive  tracking  system 
must  be  able  to  measure  and  detect  changes  in  its  performance  index  and  respond  to 
these  changes  by  performing  structural  changes.  These  structural  changes  may  involve  the 
sensing  resources  as  well  as  the  tracking  algorithms  [20] .  Chapter  4  will  describe  an  adaptive 
tracking  strategy  where  feedback  control  is  applied  on  sensing  resources  in  order  to  achieve 
given  tracking  objectives. 
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3  Discrimination  power 


Discrimination  power  is  defined  as  the  capability  to  distinguish  two  or  more  observed  ob¬ 
jects.  State  estimation  is  at  the  source  of  the  discrimination  power  problem.  If  true  states 
could  be  obtained  instead  of  estimates,  discrimination  power  would  not  be  an  issue.  The 
next  paragraphs  will  discuss  the  importance  of  discrimination  power  in  target  tracking 
applications. 

A  tracking  system  is  mainly  characterized  by  the  sensors  and  the  algorithms  that  process 
the  sensor  measurements  to  provide  state  estimation  of  targets.  Sensors  have  two  principal 
characteristics  that  come  into  play  in  the  discrimination  power  problem.  These  are  accuracy 
and  resolution,  as  illustrated  in  Figure  5. 


Figure  5:  Sensor  with  accuracy  a  (Gaussian  probability  distribution)  and  resolution  ip 

Note  that  sensor  accuracy  and  estimate  accuracy  are  defined  differently.  In  case  of  a  sensor, 
the  accuracy  is  defined  according  to  the  measurement  error  probability  distribution,  which 
is  assumed  Gaussian.  Practically,  the  value  of  accuracy  specified  in  most  sensing  systems 
represents  the  standard  deviation  of  the  distribution.  On  the  other  hand,  the  accuracy 
associated  with  an  estimation  represents  its  uncertainty  with  respect  to  the  true  value. 
Under  the  normality  assumption,  it  is  given  by  the  standard  deviation  of  the  estimate 
posterior  probability  distribution  function.  Equivalently,  accuracy  can  also  be  expressed  as 
the  standard  deviation  of  the  estimate’s  normally  distributed  random  error  with  0.0  mean. 
Note  that  accuracy  is  often  used  to  calculate  the  coverage  interval  of  estimates,  that  will  be 
discussed  in  Subsection  3. 2. 1.2.  Furthermore,  the  resolution  of  a  sensor  has  to  do  with  its 
ability  to  distinguish  between  objects  according  to  the  received  signals.  For  a  given  sensor, 
its  specified  value  of  resolution  is  equal  to  the  closest  distance  for  which  it  can  distinguish 
two  objects.  Therefore,  resolution  in  itself  represents  some  sort  of  discrimination  power.  It 
expresses  discrimination  power  at  the  signal  level. 

Moreover,  once  measurements  and  estimates  are  produced,  the  capability  to  distinguish 
the  different  observed  objects  may  still  be  an  issue.  In  such  a  case,  the  discrimination 
power  has  to  do  with  accuracy  rather  than  resolution.  Although  accuracy  depends  on 
the  signal  characteristics,  it  is  expressed  in  terms  of  the  statistical  characteristics  of  the 
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observations,  i.e.  the  probability  distribution  of  the  true  state  according  to  the  observations. 
The  discrimination  power  is  then  the  capability  to  distinguish  between  objects  according  to 
their  statistical  characteristics.  This  capability  is  directly  related  to  the  degree  of  similarity 
between  the  statistical  characteristics  of  the  objects.  The  more  the  statistics  of  observed 
objects  are  similar,  the  less  is  the  discrimination  power. 


According  to  this,  it  turns  out  that  resolution  and  accuracy  each  relate  to  a  different 
discrimination  power  level,  that  are  the  signal  level  and  the  statistical  level,  respectively. 
Figure  6  illustrates  each  of  the  two  levels  with  respect  to  accuracy  and  sensor  resolution. 


|  State  estimates 


Tracking 

system 


Discrimination  power 
accuracy  =  .  . 
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Measurements 
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Discrimination  power 
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Figure  6:  Discrimination  power  levels 

3.1  Signal  level  discrimination  power 

As  mentioned  above,  signal  level  discrimination  power  is  defined  as  the  capability  to  dis¬ 
tinguish  two  or  more  objects  according  to  the  signals  produced  by  the  sensors.  Suppose  for 
example,  a  one-dinrensional  (ID)  sensor  that  has  an  accuracy  represented  by  a  Gaussian 
distribution  with  a  standard  deviation  a.  The  sensor  also  has  a  resolution  <p,  which  is  the 
minimum  separation  interval  at  which  two  targets  can  be  separated.  Given  the  ID  sen¬ 
sor  and  two  targets  in  the  ID  space,  two  situations  may  happen:  i)  the  distance  between 
the  two  targets  is  below  the  sensor  resolution,  shown  in  Figure  7  (a),  and  ii)  the  distance 
between  the  two  targets  is  larger  than  the  sensor  resolution,  shown  in  Figure  7  (b). 

When  the  distance  between  the  two  targets  is  below  the  sensor  resolution,  the  sensor  will 
report  only  one  observed  target  since  it  cannot  distinguish  the  two  closely  spaced  targets 
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Figure  7:  Two  targets  observed  in  a  ID  space 

on  the  received  signal.  In  such  a  case,  the  discrimination  power  problem  has  to  do  entirely 
with  signal  processing  and  sensor  resolution;  thereby  only  better  sensor  resolution  could 
improve  the  discrimination  power.  Such  a  problem  will  not  be  discussed  further  in  this 
work. 


3.2  Statistical  level  discrimination  power 

When  the  distance  between  two  targets  is  larger  than  the  sensor  resolution  (Figure  7  (b)), 
the  sensor  will  usually  report  two  different  target  observations.  In  that  case,  the  sensor 
has  enough  discrimination  power  at  the  signal  level.  However,  even  though  a  different 
observation  is  produced  for  each  of  the  two  targets,  there  may  be  a  discrimination  power 
problem  related  to  the  estimation  process,  i.e.  the  accuracy  of  the  estimates.  As  defined  in 
Chapter  6,  the  accuracy  of  an  estimation  depends  on  the  related  probability  distribution. 
In  case  of  two  or  more  estimations,  the  closer  the  related  probability  distributions  are,  the 
lower  will  be  the  discrimination  power.  This  is  shown  in  Figure  8. 


(a)  High  discrimination  power  -  There  is  almost  (b)  Low  discrimination  power  -  There  is  a  large 
no  overlap  between  the  probability  distributions  overlap  between  the  probability  distributions  of 
of  the  estimated  states.  the  estimated  states. 


Figure  8:  Discrimination  power  for  two  state  estimates 
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Within  the  tracking  process,  filtering,  the  track  updates  ( i.e .,  state  updates  and  time  up¬ 
dates)  presented  in  Subsection  2.2.1  influence  the  accuracy  of  the  state  estimation,  which 
has  a  direct  impact  on  the  discrimination  power  over  time.  In  Chapter  4,  a  method  for 
controlling  the  discrimination  power  is  presented,  where  track  update  periods  are  adjusted 
adaptively. 

3.2.1  Discrimination  power  metrics 

Any  metric  derived  to  measure  the  discrimination  power  level  will  have  to  be  based,  in  one 
way  or  another,  on  the  probability  distributions  of  objects  in  space  and  time.  Therefore, 
before  deriving  such  metrics,  it  is  important  to  recall  some  basic  concepts  around  the  notion 
of  probability. 

3.2. 1.1  Probability 

There  exist  different  interpretations  of  what  probability  is  about.  Generally,  statisticians 
divide  into  two  groups  in  their  approach  for  interpreting  probability,  the  Bayesians  and  the 
frequentists  [23,  24,  25,  26,  27]2. 

With  the  Bayesian  approach,  probability  is  interpreted  as  a  degree  of  belief  that  a  parameter 
will  have  a  given  value.  The  degree  of  belief  is  conditional  on  all  relevant  information  about 
the  parameter.  Practically,  this  assumes  a  prior  distribution  for  the  parameter  to  estimate. 
Bayesian  statistics  are  the  basis  for  target  tracking  algorithms  such  as  Kalman  filtering. 

Often  said  to  be  closer  to  scientific  reasoning  [27],  the  frequentist  approach,  also  called 
conventional  or  frequency-based  approach,  interprets  probability  as  the  relative  frequency 
of  something  happening.  Frequentists  consider  parameter  estimation  as  the  result  of  a 
repetition  of  experiments,  where  an  estimated  parameter  is  a  random  variable  that  can 
take  on  different  values  when  the  experiment  is  repeated.  On  the  other  hand,  Bayesian 
statistics  do  not  consider  estimation  as  a  repetition  of  experiments,  where  an  estimated 
parameter  is  seen  as  an  unknown  variable  instead  of  a  random  variable. 

Depending  on  the  approach,  Bayesian  or  frequentist,  measuring  the  discrimination  power 
may  have  different  meanings.  We  present  two  ways  of  measuring  the  discrimination  power, 
coverage  intervals3  and  Mahalanobis  distance.  The  former  refers  to  the  Bayesian  approach, 
while  the  latter  is  a  frequentist  method.  The  measures  consider  two  estimates  at  a  time  for 
calculating  the  discrimination  power. 

2Note  that  R.A.  Fisher  presented  an  approach  for  interpreting  statistics  that  is  also  popular  among  the 
scientific  community  [28,  29].  In  short,  his  approach  consists  of  a  series  of  theories  that  offer  compromises 
between  the  Bayesian  and  frequentist. 

3Note  that  coverage  intervals  are  to  the  Bayesian  approach  what  confidence  intervals  are  to  the  frequentist 
approach. 
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3.2. 1.2  Coverage  intervals 


A  Coverage  Interval  (Cl)  is  an  interval  on  the  probability  distribution  function  associated 
with  the  estimated  state.  It  has  a  different  interpretation  from  the  confidence  intervals  in 
the  frequentist  approach  [24,  25]. 

Suppose  a  state  x  is  estimated  from  some  measurements  Zk  =  [z\,Z2,  —  ,Xk\  and  that  the 
resulting  probability  distribution  is  N(x,a).  A  coverage  interval  is  an  interval  obtained 
from  the  distribution  N(x,  a)  that  expresses  how  much  the  estimated  value  x  can  at  most 
deviate  from  the  unknown  true  value  x  [30].  The  interval  has  an  associated  probability  a 
that  it  includes  the  true  value,  also  called  degree  of  belief.  Generally,  the  recommended 
degree  of  belief  is  a  =  95%,  which  is  illustrated  in  Figure  9  and  corresponds  to  2  standard 
deviations  (2a)  of  the  estimated  value  for  a  Gaussian  distribution  [31]. 


Figure  9:  95%  coverage  interval  (Cl)  of  an  estimate  x 

Under  the  Bayesian  approach,  the  interval  shown  in  Figure  9  means  that  chances  are  95% 
that  the  true  value  of  x  be  included  in  the  interval4. 

Considering  two  different  estimates  x\  and  X2  from  group  of  measurements  Zk  and  Zk, 
respectively,  the  discrimination  power  will  be  determined  according  to  their  respective  CIs, 
CI\  and  CI2 ,  for  a  degree  of  belief  a  with 

CI\  =  [Ci,  C?]  (8) 

and 

Ch  =  [ClC%]  (9) 

and  where  l  stands  for  “lower  limit”  and  u  stands  for  “upper  limit”.  Thus,  when  the 
coverage  intervals  of  two  estimates  do  not  overlap,  that  is 

[Cl,  Ci]  C  [Cl2,C%]  =  0,  (10) 

4With  the  frequentist  approach,  the  interpretation  is  that  the  calculated  confidence  interval  has  a  95% 
probability  of  containing  the  true  states  x.  In  other  words,  if  confidence  intervals  were  calculated  repeatedly 
with  new  samples  (or  measurements),  95%  of  them  would  contain  the  true  state  x. 
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it  is  assumed  that  there  is  a  statistically  significant  difference  between  the  estimates  x±  and 
X‘2-  It  is  then  understood  that  the  two  estimated  objects  can  be  clearly  distinguished  and 
that  the  discrimination  power  is  fair.  This  is  shown  on  Figure  10  where  Cf  =  C2. 


Figure  10:  Confidence  intervals  and  discrimination  power  -  CI\  and  CI-2  do  not  overlap, 
the  discrimination  power  is  fair 


On  the  contrary,  when  the  coverage  intervals  of  two  estimates  overlap,  that  is 

[c[,cf]  n  [cl2,c%]  +  0,  (li) 

the  conclusion  is  that  there  is  not  a  statistically  significant  difference  between  the  estimates 
x\  and  S''}. ■  Therefore,  the  two  estimated  objects  cannot  be  clearly  distinguished  and  the 
discrimination  power  is  too  low.  This  is  shown  on  Figure  11,  where  Cf  >  C2.  Note  that 
under  the  frequentist  approach,  it  is  not  necessarily  true  that  there  is  not  a  statistically 
significant  difference  between  both  estimates  when  the  confidence  intervals  do  overlap  [32] . 
Hence  this  method  does  not  fit  with  the  frequentist  viewpoint. 


Figure  11:  Confidence  intervals  and  discrimination  power  -  CI\  and  CI2  overlap,  the  dis¬ 
crimination  power  is  too  low. 


Although  coverage  intervals  are  easy  for  the  ID  case,  the  method  is  more  complex  for 
the  2D  (and  more)  case.  In  2D  for  example,  applying  the  coverage  intervals  concept  would 
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require  calculating  intersections  of  probability  distributions  contours  (ellipsoids  for  Gaussian 
distributions).  Other  methods  should  then  be  considered,  such  as  the  Mahalanobis  distance 
measure. 


3.2. 1.3  Mahalanobis  distance 


The  Mahalanobis  distance  is  a  statistics  that  computes  a  difference  between  two  estimates 
by  taking  into  account  their  variances  [33].  The  Mahalanobis  distance  d  between  two 
estimates  X{  and  Xj  is  defined  as. 


T 

-1 

- 

d(xi  ?  Xj )  — 

Xi  —  Xj 

Pi  +  Pj 

Xi  —  Xj 

(12) 


where  Pj  and  Pj  are  the  covariance  matrices  for  the  state  estimates  x j  and  Xj,  respectively. 
Because  of  the  Gaussian  assumption,  d(xi,Xj )  is  a  \2  random  variable  with  degrees  of 
freedom,  where  is  the  number  of  dimensions  for  the  state  estimates  Xi  and  Xj. 


For  N  estimates,  the  distances  for  each  possible  pair  of  estimates  are  represented  in  the 
N  x  N  dissimilarity  matrix 


V  = 


d(xi,Xj) 


(13) 


The  Mahalanobis  distance  can  be  used  as  a  measure  of  the  similarity  between  tracks.  Fig¬ 
ure  12  illustrates  pairs  of  tracks’  estimates  with  different  levels  of  similarity,  represented  by 
the  statistical  distance  d. 


Figure  12:  Statistical  distance  d  as  a  track  dissimilarity  metric 


The  ellipsoids  represent  contours  of  constant  probability  for  the  two-dimension  and  normally 
distributed  position  estimates  of  the  tracks.  As  can  be  seen  in  Figure  12,  the  farther 
the  position  estimates  x\  and  X2  are  from  each  other  in  the  space,  the  higher  are  the 
Mahalanobis  distance  and  track  dissimilarity  (d±  >  c^)-  Also,  the  smaller  the  estimation 
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error  covariance  matrices  P\  and  P2  are,  the  higher  are  the  Mahalanobis  distance  and  track 
dissimilarity  is  (g?2  >  c/3). 

3.2.1 .4  x2  test  for  assessing  discrimination  power 

Furthermore,  given  the  Mahalanobis  distance,  how  would  one  determine  whether  the  dis¬ 
crimination  power  is  fair  or  not?  To  answer  this  question,  let  us  consider  the  following 
hypotheses  for  the  two  estimates  Xjt  and  xy. 

H0.  x i  and  Xj  refer  to  the  same  object.  This  hypothesis  suggests  that  the  discrimination 
power  is  not  high  enough. 

Hr  This  is  the  alternative  hypothesis  that  xt  and  Xj  refer  to  different  objects.  This 
hypothesis  suggests  that  the  discrimination  power  is  sufficiently  high. 

Statistically,  testing  hypothesis  Hq  against  H\  comes  down  to  determine  whether  estimates 
Xi  and  Xj  come  from  the  same  distribution  or  not.  The  Mahalanobis  distance  evaluated  in 
Equation  12  is  used  to  determine  whether  Hq  should  be  rejected  (and  replaced  with  H\ ) 
or  not.  The  criterion  for  deciding  between  Hq  or  H\  will  be: 

if  d(xi,Xj)  <  X2(a,ndy,  Accept  H0  (14) 

if  d(xi,sbj )  >  X(Qnd)>  Reject  Hq  and  accept  H\  (15) 

where  xfa  „d)  is  the  chi-square  percent  point  function  with  nd  degrees  of  freedom  and 
a  significance  level  of  a.  The  significance  level  a  represents  the  confidence  in  hypothesis 
Hq  [34,  35].  That  is,  if  Hq  is  true,  there  should  be  a  1  -a  probability  that  d(xi ,  Xj)  <  xfa  ndy 
Thus,  the  lower  the  significance  level  a,  the  higher  the  confidence  should  be  about  accepting 
or  rejecting  Hq.  In  terms  of  discrimination  power,  this  significance  level  a  represents  a 
measure  of  the  degree  of  discrimination  power.  As  a  decreases,  the  corresponding  expected 
discrimination  power  level  increases. 

Based  on  the  definition  of  discrimination  power  presented  in  the  current  chapter,  the  next 
chapter  proposes  an  approach  to  control,  adaptively,  the  discrimination  power  over  a  certain 
volume  of  interest. 
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4  Discrimination  power  control 


This  chapter  discusses  the  problem  of  controlling  the  discrimination  power,  at  the  statistical 
level,  in  a  tracking  system.  The  goal  is  to  maximize  the  discrimination  power  of  the  tracks 
in  a  certain  volume  of  interest,  such  that  all  tracks  can  be  clearly  distinguished.  This  work 
concentrates  on  controlling  the  discrimination  power  by  adjusting  track  update  periods, 
that  represent  time  intervals  between  sensor  measurements. 

First,  a  simple  algorithm  for  adaptively  updating  tracks  is  presented  in  Section  4.1.  Then, 
a  clustering-based  control  method,  that  is  better  suited  for  target  tracking  systems,  is 
presented  in  Section  4.2. 


4.1  Simple  track  update  algorithm 

A  solution  to  control  the  track  dissimilarity  (and  the  related  discrimination  power)  consists 
in  adjusting  adaptively  the  time  intervals  between  updates  for  the  different  tracks.  The 
time  interval  between  measure  updates  will  also  be  referred  to  as  the  update  period.  It 
has  a  significant  influence  on  the  variation  in  time  of  the  track  accuracy  ( i.e .,  covariance 
matrix) .  Shorter  periods  between  updates  should  yield  higher  track  accuracy  and  therefore 
higher  dissimilarity  between  tracks. 

Therefore,  the  goal  of  an  adaptive  tracking  system  is  to  choose  appropriate  update  periods 
for  each  track  in  order  to  keep  a  certain  level  of  discrimination  level  within  the  volume  of 
interest. 


A  straightforward  discrimination  power  control  solution  is  to  have  each  track’s  own  up¬ 
date  period  h  be  a  function  of  the  distance  with  its  nearest  neighbor,  according  to  the 
Mahalanobis  distance.  A  track  i  would  have  an  update  period  hi  determined  such  that 


hi 


c 


min  d(xi,xk)  , 
k 


k  el,  ...,N 


(16) 


where  d(xi,Xj )  is  as  defined  by  Equation  12  and  £  is  a  function  that  expresses  the  relation 
between  h  and  d. 


Resolving  Equation  16  would  be  complicated  practically,  requiring  a  lot  of  computations. 
Furthermore,  the  resulting  strategy  could  hardly  match  with  the  scheduling  of  sensors. 


4.2  Clustering-based  control  method 

The  clustering-based  method  aims  at  separating  the  volume  of  interest  into  regions  having 
different  discrimination  power  levels.  The  tracks  in  each  region  are  then  updated  at  some 


DRDC  Valcartier  TR  2008-016 


17 


frequencies  that  will  help  enhance  the  discrimination  power  in  regions  where  the  latter  is 
low.  The  method  relies  on  an  adaptive  control  structure  suited  for  the  specific  problem  of 
target  tracking  and  depicted  in  Figure  13.  It  is  a  practical  implementation  of  the  adaptive 
data  fusion  concepts  presented  in  Section  2.1  and  illustrated  in  Figure  4. 


Mission 

Objectives 


Figure  13:  Clustering/discrimination  power-based  Two-Level  Control  Loop 


The  control  structure  is  made  of  a  two-level  cascade  control  loop.  Its  two  loops  operate  at 
two  different  time-scales,  with  the  following  objectives 

•  Outer-loop  —  The  characterization  of  the  whole  volume  of  interest  based  on  the 
concept  of  discrimination  power,  where  distinct  regions  are  created  and  managed  dy¬ 
namically  based  on  the  targets  spatial  distribution  and  by  using  clustering  techniques. 
To  ensure  a  certain  level  of  cluster  persistence,  this  loop  operates  at  low  frequency, 
as  defined  by  the  basic  period  hc  >  hs  in  Figure  13.  Thus  sensor  scheduling  can  be 
made  by  considering  several  sampling  intervals. 

•  Inner-loop  —  The  adaptive  control  of  discrimination  power  over  the  volume  of  inter¬ 
est,  by  the  appropriate  selection  of  the  track  update  frequency  and  sensor  scheduling 
within  each  distinct  region.  This  adaptation  loop  operates  at  a  faster  time  scale  than 
the  outer  loop  (hs  <  hc).  Note  that  the  minimal  value  for  hs  is  a  characteristic  of  the 
sensors. 

4.2.1  Hierarchical  clustering  of  tracks 

When  regions  of  the  VOI  are  to  be  considered  for  controlling  sensor  scans,  separating  tracks 
into  clusters  can  provide  a  way  to  determine  the  desired  update  strategy.  Clusters  of  targets 
represent  different  regions  of  the  space  with  different  discrimination  power  levels. 

Agglomerative  hierarchical  clustering  [33]  is  suitable  for  controlling  the  discrimination  power 
using  the  region-based  approach.  With  agglomerative  hierarchical  clustering,  a  binary  tree 
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is  constructed  based  on  a  linkage  measure,  where  each  leaf  represents  the  predicted  state 
estimate  of  a  track  (Figure  14).  The  linkage  measure  is  determined  and  related  to  the  type 
of  hierarchical  clustering  algorithm. 


9 

9 


Figure  14:  Distance  tree  and  hierarchical  formation  of  clusters 

A  hierarchical  cluster  tree  is  created  based  on  the  dissimilarity  matrix  V.  The  process  of 
hierarchical  clustering  is  as  follows: 

1.  Assign  each  track  to  a  cluster  to  have  N  clusters. 

2.  Merge  the  closest  pair  of  clusters,  where  cluster  Ci  and  Cj  are  merged  together  to 
result  in  N  —  1  clusters. 

3.  Compute  distances  between  the  new  cluster  and  each  of  the  old  clusters.  The  distance 
computation  depends  on  the  type  of  clustering  used.  The  following  gives,  for  two 
clusters  Ci  and  Cj,  the  resulting  distances  dmin(Ci,  Cj)  and  dmax{Ci,  Cj),  when  single- 
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linkage  and  complete-linkage  are  used,  respectively 


dmin  (  Ci ,  Cj  ) 

=  min  d(x,  x  ) 

(17) 

x£Ci,x'£Cj 

dmax  (  Ci  i  Cj  ) 

=  max  d(x,x') 
xeCi,x'eCj 

(18) 

4.  Repeat  Steps  2  and  3  until  all  tracks  are  clustered  into  a  single  cluster  of  size  N. 

The  single-linkage  algorithm  produces  a  binary  tree  where  the  highest  nodes  in  the  hierarchy 
have  the  largest  distance  values.  A  threshold  L  applied  on  the  distances  allows  to  group 
tracks  into  regions.  The  value  of  the  threshold  determines  how  the  tree  is  separated.  All 
nodes  that  are  below  the  threshold  will  have  their  corresponding  tracks  regrouped  into  the 
same  cluster.  Each  cluster  of  tracks  occupies  a  particular  region  of  the  space.  The  value 
of  the  threshold  should  depend  on  the  problem  at  hand,  that  is  the  number  of  targets 
and  their  properties,  the  number  of  sensors  and  their  properties  and  most  importantly  the 
applications  that  exploits  the  results  of  the  target  tracking  process5. 

4.2.2  Sensor  scheduling  heuristics 

Once  the  clusters  of  targets  are  determined,  the  scheduling  of  the  sensors  is  required.  Here, 
a  single  sensor  is  considered.  It  is  assumed  that  the  sensor  can  be  allocated  to  one  cluster 
at  a  time.  It  is  also  assumed  that  the  sensor  can  be  directed  instantaneously,  without  extra 
cost. 


Let  Rn  be  the  region  occupied  by  cluster  Cn,  for  which  the  discrimination  power  dn  and 
the  track  similarity  sn  are  defined  as  follows 


sn  =  min d(xi(k),Xj(k)),  &i,Xj  6  Rn  (19) 

dn  =  s-1  (20) 

The  two  metrics  dn  and  sn  will  be  used  to  define  the  optimal  track  update  frequency  and 
the  scheduling  strategy.  The  update  cycle  length  of  the  clusters  is  defined  as 


hc 


hs 

I< 

hs£  E  7?n 


n= 1 
K 


hstJ2  L— 

•  m  l  r 


d 


-i 


d 


i=i  mmi6{lr  ^q 

^El- — - — . 

^  mmie{lr..)K}  Si 


-l- 


5e.g.,  surveillance,  target  engagement,  or  other. 


(21) 

(22) 

(23) 

(24) 
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where  K  is  the  number  of  clusters  and  £  is  an  integer  defined  arbitrarily.  The  integers  r]n 
define  the  number  of  updates  that  will  be  performed  for  each  region  within  a  single  update 
cycle  i p.  The  latter  is  defined  such  that  the  tracks  within  the  cluster  that  has  the  highest 
discrimination  power  ( i.e smallest  sn)  will  be  updated  only  once,  while  the  tracks  within 
the  other  clusters  will  be  more  than  once. 

It  is  assumed  that  the  sensor  has  to  spend  a  minimum  time  ts  over  a  region  to  report 
a  contact.  This  defines  a  maximum  report  frequency  fs.  The  update  frequency  for  each 
cluster,  as  function  of  the  sensor  frequency  fs,  is  defined  as  follows, 

fn  =  C V“  Vis  (25) 

For  regions  that  require  more  than  one  update,  different  update  scheduling  strategies  may  be 
possible.  Three  different  scheduling  algorithms  are  considered  and  compared  in  the  sequel. 
These  are:  time-slice,  minimum  intra-cluster  discrimination  power  and  round-robin. 

4.2.2. 1  Time-slice 

This  strategy  is  described  by  Algorithm  1,  where  the  objective  is  to  separate  (in  time)  the 
updates  of  the  same  region  as  much  as  possible  in  order  to  maintain  a  best  quality  of  track 
with  the  same  number  of  updates,  as  illustrated  in  Figure  15.  The  schedule  period  length 
corresponds  to  the  cluster  update  cycle  length  hc. 

1:  procedure  TimeSlice 

2:  S  <—  schedule  vector  of  size  V’  >  Update  cycle  length 

3:  j  <-  k  |  rjk  =  max/e{li...iA'}  rji 

4:  i  <—  1 

5:  while  maxigp  Aq  r/i  >  0  do 

6: 

7:  Vk  <-  r)k  ~  1 

8:  if  rjk  >  1  then 

9:  S(lp  —  i  +  1)  <—  j 

10:  rjk  rik  -  1 

11:  end  if 

12:  i<—i+  1 

13:  j  <-  k  |  rik  =  nrax;e|lj  A-|  //; 

14:  end  while 

15:  end  procedure 

Algorithm  1.  Time-Slice  Scheduling  Strategy 
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No  Update 


Figure  15:  Track  Quality  as  function  of  scheduling  strategy  for  the  same  number  of  updates 

4.2.2.2  Minimum  Intra-Cluster  Discrimination  Power  (MICDP) 

With  the  minimum  intra-cluster  discrimination  power  scheduling  algorithm,  the  next  region 
Rm  to  be  updated  is  the  one  with  the  minimum  discrimination  power: 

m  I  dm  =  min  dn  (26) 

1  ne{i 

where  the  clusters  are  being  updated  after  each  sensor  update,  that  is 

hc  =  K  (27) 


4.2.2.3  Round-robin 

The  round-robin  scheduling  strategy  updates  all  the  clusters  equally.  With  such  a  strategy, 
a  cluster  that  has  just  been  updated  will  not  be  updated  again  until  all  other  clusters  have 
been  updated.  The  update  cycle  length  of  the  clusters  is  the  same  as  for  the  time-slice 
strategy  defined  in  Equation  24.  The  round-robin  strategy  represents  the  open-loop  and 
static  approach  for  controlling  the  sensing  operations. 
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5  Results  and  Discussion 


Several  scenarios,  in  which  a  single  (phased  array-like)  sensor  [36]  has  to  track  several 
targets,  were  simulated.  For  each  of  the  scenarios,  different  scheduling  strategies  were 
compared  on  the  base  of  different  performance  measures.  The  following  sections  describe 
the  two-dinrension  phased  array  radar  model  and  the  performance  measures  that  are  used 
in  the  simulation  in  order  to  test  and  compare  the  different  scheduling  strategies. 

5.1  Sensor  Model 

The  sensor  model  assumes  an  active  sensor  that  has  the  capability  to  switch  the  direc¬ 
tion  of  its  beam  instantaneously  without  inertia.  Precisely,  the  sensor  has  the  following 
characteristics: 

1.  The  sensor  measures  range  and  azimuth. 

2.  The  beam  width  in  elevation  is  considered  large  enough  to  cover  the  entire  elevation 
span.  Therefore,  this  is  a  2D  sensor. 

3.  The  beam  direction  is  controllable  in  azimuth. 

4.  Measures  are  made  instantaneously.  Hence,  the  search  time  for  a  track  update  is  zero. 

5.  Detection  probability  is  100%. 

6.  False  alarm  probability  is  0%. 

7.  In  some  cases,  it  is  assumed,  for  simplicity,  that  more  than  one  measure  on  different 
targets  can  be  made  simultaneously.  That  is,  all  targets  pertaining  to  a  same  cluster 
are  updated  instantaneously. 

8.  Standard  deviations  of  range  and  azimuth  measurements  are  assumed  constant.  They 
are  55.6m  and  7t/90  rad  (or  0.5°)  in  range  and  azimuth,  respectively. 

5.2  Performance  measures 

To  evaluate  the  performance  of  the  proposed  solution,  several  metrics  have  been  defined. 
They  are  all  function  of  the  uncertainty,  as  expressed  by  the  error  covariance  matrix  of  each 
state  estimate.  A  sub-set  of  these  metrics  is  described  below. 
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5.2.1  Discrimination  Conflict 


We  assume  that  there  is  a  Discrimination  Conflict  (DC)  between  two  different  state  esti¬ 
mates  x,  and  Xj,  at  time  k,  if  their  Mahalanobis  distance  [33]  is  lower  than  (or  equal  to)  a 
certain  threshold  A^,  that  is 


DCij(k) 


1  if  (xi(k),Xj(k))  <  Xd 
0  otherwise 


(28) 


Thus,  there  is  a  discrimination  conflict  if  the  overlap  between  the  spatial  probability  distri¬ 
butions  of  two  targets  exceeds  some  threshold  value. 


Discrimination  conflicts  should  stay  as  low  as  possible  to  avoid  problems  with  contact/track 
correlation,  target  identification  and  classification.  The  total  number  of  discrimination 
conflicts  of  a  tracking  system  over  a  period  of  time  [1, ...,  M]  is  the  sum  of  the  conflicts  for 
each  discrete  time  steps  k: 

M 

DCit,,,tM  =  EE  DCij(k),  i,jel,...,N  (29) 

k= 1  j>i 


5.2.2  Average  track-to-track  Mahalanobis  distance 

The  average  Mahalanobis  distance  for  any  pairs  of  tracks  within  a  given  cluster,  at  a  given 
time  k,  is  defined  as 

AM(k)=  v  y  i,jel,...,N  (30) 

This  metrics  is  a  characteristic  of  the  clusters.  The  higher  is  the  average  track-to-track 
Mahalanobis  distance,  the  better  is  the  discrimination  power  within  the  clusters. 

5.2.3  Average  error  covariance  matrix  determinant 

The  average  error  covariance  matrix  determinant  at  a  time  k  is  defined  as 

1  N 

ADET(k)  = -Y,\pi(k)\  (31) 

i= 1 

This  measure  provides  a  means  to  express  the  overall  (for  the  all  volume  of  interest)  uncer¬ 
tainty  of  the  tracks  at  a  specific  time. 


5.3  Scenarios 

Three  scenarios  (A,B,C)  are  simulated  and  the  three  metrics  defined  above  are  used  to 
evaluate  the  performance  of  the  proposed  solution.  Scenarios  A  and  B  feature  pre-scripted 
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trajectories  of  four  (4)  and  seven  (7)  targets,  respectively,  while  Scenario  C  uses  fifteen  (15) 
randomly  distributed  targets. 

The  goal  of  the  control  loop  for  the  following  scenarios  is  to  maintain  the  appropriate 
dissimilarity  level  between  the  tracks’  latest  estimates.  The  control  strategy  should  make 
the  sensor  spend  more  time  over  those  regions  where  targets  have  low  dissimilarity. 

The  Euclidian  distance  threshold  for  creating  clusters  is 

L  =  400m  (32) 

The  Mahalanobis  distance  threshold  for  calculating  the  total  number  of  discrimination 
conflicts  is 

Xd  =  50  nr  (33) 

where  Xd  corresponds  to  the  xfQ  n({)  vahie  discussed  in  Sub-section  3. 2. 1.4.  The  update  cycle 
length  of  the  clusters  hc  and  the  track  update  period  hs  are  1.0  s  and  0.1  s  respectively. 


The  targets  move  with  an  assumed  constant  velocity  model.  The  process  noise  v k  defined 
in  Equation  1  is  white  with  zero-mean  and  has  a  covariance 


E(vkv'k)  = 


(34) 


with 

q  =  625 m2/s3 

and  where 

7i2/2  0 

0  h2/ 2 

h  0 

0  h 


(35) 


(36) 


for  a  time  interval  h 
h  is 


t'k+i  —  ffc-  The  state  transition  matrix  Ff.  for  the  same  time  interval 


Fk 


1  0  h  O' 
0  10  h 
0  0  10 
0  0  0  1 


(37) 


The  phased  array  radar  measurements  are  converted  to  Cartesian  coordinates  using  the 
conventional  coordinate  transformation  [3,  37]. 
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5.3.1  Scenario  A 


To  illustrate  the  benefits  of  the  discrimination  power  control  strategy,  a  simple  scenario 
with  only  four  (4)  targets  has  been  simulated,  where  two  of  the  four  targets  cross  each 
other  at  a  given  time.  Figure  16  illustrates  the  trajectories  of  the  4  targets.  Targets  1  and 
2  are  set  to  cross  at  about  5  s,  while  targets  3  and  4  are  set  far  from  each  other.  Following 
this  scenario,  it  is  expected  that  the  tracking  system  will  have  more  difficulty  discriminating 
Target  1  and  Target  2  during  their  intersection  period. 


I 

3 

xIO4 


Figure  16:  Trajectories  of  the  4  targets  in  Scenario  A  (‘X’  is  the  final  position) 

Figure  17  shows  the  dissimilarity  of  the  tracks  corresponding  to  Target  1  and  Target  2. 
The  dissimilarity  is  expressed  by  the  Mahalanobis  distance6.  The  Mahalanobis  distance 
is  shown  for  the  first  8  seconds  of  the  simulation  and  the  different  scheduling  strategies 
presented  in  Sub-sections  4. 2. 2.1,  4. 2. 2. 2  and  4. 2. 2. 3.  For  the  period  of  time  where  targets 
1  and  2  are  close  to  each  other,  that  is  between  about  3  s  and  7  s  from  the  start  of  the 
simulation,  the  Mahalanobis  distance  gets  smaller  and  reaches  a  minimum  at  5  s,  that  is 
when  targets  1  and  2  are  crossing.  As  illustrated  on  Figure  17,  the  time-slice  scheduling 
strategy  is  the  one  that  was  able  to  obtain  the  best  discrimination  power  by  maintaining 
the  highest  Mahalanobis  distance,  compared  to  the  minimum  intra-cluster  discrimination 
power  (MICDP)  and  round-robin  scheduling  strategies. 

5.3.2  Scenario  B 

A  scenario  with  seven  (7)  targets  has  been  simulated  and  the  results  are  presented  in  this 
section.  Within  the  considered  configuration,  certain  targets  will  come  too  close,  at  given 
time  instants,  to  be  clearly  distinguishable  from  each  other,  e.g.  Target  1  and  Target  2  at  5 
s,  and  Target  1  and  Target  3  at  14  s,  as  shown  on  Table  1  where  the  closest  pairs  of  targets 
are  given  in  terms  of  different  simulation  times. 

Figures  18  and  19  present  and  compare  the  resulting  track  quality  based  on  two  different 
track  update  strategies,  for  a  tracking  duration  of  20  s.  Dashed  ellipses  give  initial  position 

6Note  that  any  well-defined  distance  on  Rn  may  be  used  as  a  proximity  metric.  The  presented  results  and 
the  underlying  development  are  based  on  the  Mahalanobis  distance  since  it  considers  the  error  covariances. 
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- Round  robin  scheduling 

. Time-slice  scheduling 

- MICDP  scheduling 

- Ground  Truth 
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Figure  17:  Mahalanobis  distance  between  tracks  of  targets  1  and  2.  The  ground  truth 
distance  for  targets  1  and  2  is  also  illustrated. 


Time  (s) 

Closest  pair 

Distance  (nr) 

0 

1  -  2 

1300 

5 

1  -  2 

10 

10 

1  -  3 

1400 

14 

1  -  3 

5 

20 

1  -  4 

1100 

Table  1:  Closest  pairs  of  targets  in  scenario  B. 


and  uncertainty  for  each  target.  Plain  ellipses  are  represented  only  to  show  final  position 
and  uncertainty  for  each  target,  and  also  where  discrimination  problems  are  expected,  e.g., 
intersections  (Target  1  [red],  Target  3  [green])  and  (Target  1  [red],  Target  2  [blue]). 

In  Figure  18,  a  round-robin  update  strategy  was  used,  where  the  different  clusters  of  tracks 
are  being  allotted  the  same  attention  without  consideration  of  any  additional  information. 
On  the  contrary,  Figure  19  presents  the  time-slice  strategy,  the  objective  of  which  is  the 
control  of  the  discrimination  power  within  the  different  clusters.  More  attention  is  paid 
to  critical  areas,  i.e.,  where  discrimination  power  tends  towards  zero.  Clusters  are  created 
dynamically  based  on  spatial  proximity,  and  the  update  strategy  is  defined  based  on  the 
track  dissimilarity  in  each  cluster. 

Figure  19  shows  clearly  the  superiority  of  the  time-slice  discrimination  power  control-based 
scheduling  approach  over  the  static  round-robin  scheduling  strategy.  Discrimination  power 
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is  adaptively  improved  where  required,  e.g.,  over  intersections  (Target  1  [red],  Target  3 
[green])  and  (Target  1  [red],  Target  2  [blue]),  by  increasing  the  update  rate  over  the  clusters 
created  by  the  targets  proximity.  Targets  that  do  not  need  high  discrimination  power  will 
belong  to  distinct  clusters  that  will  be  updated  less  frequently,  e.g.,  Target  6  [cyan]  and 
Target  5  [yellow].  For  the  least  frequently  updated  tracks  by  the  adaptive  strategy  (i.e., 
Target  3  [cyan]  and  Target  4  [yellow]),  the  obtained  results  show  1.32  km 2  vs.  1.0  km2  with 
static  for  Target  6  [cyan]  and  1.11  km 2  vs.  .86  km 2  with  static  for  Target  5  [yellow]. 


[km] 


Figure  18:  Track  quality  yielded  by  round-robin  scheduling  strategy 


Note  that,  in  Figures  18  and  19,  the  ellipses  that  are  not  centered  on  the  last  update 
correspond  to  tracks  to  which  sensor  time  was  not  allocated  during  the  last  part  of  the 
scenario.  For  these  tracks,  only  prediction  (dead  reckoning)  is  used. 

5.3.3  Scenario  C 

To  further  evaluate  the  performance  of  the  proposed  strategy,  a  second  scenario  with  fifteen 
(15)  randomly  distributed  targets  is  used.  For  both  axes  X  and  Y  in  the  Cartesian  plane,  the 
target  initial  positions  are  bounded  between  [1000/cm,  3000/cm].  The  target  initial  velocities 
are  determined  randomly  in  the  range  [0  m/s, 250  m/s]  over  axis  X  and  are  set  to  0  over  axis 
Y.  These  give  the  initial  velocity  values  that  may  change  during  the  simulation  because  of 
a  deliberately  added  noise  (see  Figure  20).  The  simulation  duration  is  15  s. 
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Figure  19:  Track  quality  yielded  by  the  time-slice  discrimination  power  control-based 
scheduling  strategy 

An  example  of  the  fifteen  target  trajectories  is  illustrated  in  Figure  20.  Time-slice,  minimum 
intra-cluster  discrimination  power  (MICDP)  and  round-robin  scheduling  strategies  were 
tested.  Figures  21  to  24  show  the  resulting  performance  measures  according  to  100  Monte 
Carlo  runs  for  each  of  the  three  scheduling  strategies. 

The  number  of  clusters,  averaged  over  the  100  Monte  Carlo  runs,  is  given  in  Figure  21. 
Thus,  on  average  the  number  of  clusters  is  initially  around  7  and  increases  slightly  with 
time.  Such  an  increase  is  caused  by  target  dispersion  as  time  goes  on. 

The  average  track-to-track  Mahalanobis  distance  (AM)  is  given  in  Figure  22.  It  shows 
that  the  highest  AM  is  obtained  with  the  time-slice  scheduling  strategy.  The  latter  allows 
maintaining  the  highest  discrimination  power  over  the  volume  of  interest.  The  MICDP 
strategy  performs  worse  than  the  round-robin  method  and  yields  the  lowest  AM . 

Moreover,  time-slice  scheduling  provides  also  the  best  result  when  it  comes  to  minimizing 
the  total  number  of  discrimination  conflicts,  as  shown  in  Figure  23.  The  latter  illustrates  the 
distribution  of  the  total  number  of  discrimination  conflicts  (PC^ ...,15])  over  the  100  Monte 
Carlo  runs.  Discrimination  conflicts  occur  much  more  often  with  the  MICDP  scheduling 
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Figure  20:  Sample  of  trajectories  of  the  fifteen  targets  in  Scenario  C  (‘X’  is  the  final  position) 

strategy,  whose  average  total  number  of  similarity  conflicts  is  2017,  compared  to  1400  and 
623  for  the  round-robin  and  time-slice  scheduling  methods,  respectively.  In  Figure  23,  the 
lower  and  upper  lines  of  the  box  are  the  25th  and  75th  percentiles  of  the  samples.  The  line 
in  the  middle  of  the  box  is  the  sample  median.  The  lines  extending  above  and  below  the 
box  show  the  extent  of  the  rest  of  the  samples. 

Figure  24  shows  that  maintaining  high  discrimination  power  does  not  necessarily  mean 
maintaining  an  overall  low  individual  track  uncertainty.  It  represents  the  time  evolution  of 
the  average  determinant  of  the  error  covariance  matrix  ( ADET ).  It  is  clearly  shown  that 
though  time-slice  is  the  best  scheduling  strategy  for  maintaining  high  discrimination  power, 
it  is  not  necessarily  the  method  that  maintains  the  lowest  individual  track  uncertainties 
overall.  For  example,  between  1  s  and  4  s  on  Figure  24,  round-robin  scheduling  provides  a 
lower  ADET  than  time-slice  scheduling  over  the  100  Monte  Carlo  runs. 

Although  both  time-slice  and  MICDP  aim  at  maximizing  the  discrimination  power,  the 
results  have  demonstrated  that  time-slice  scheduling  performs  better.  The  latter  focuses 
the  sensing  effort  where  and  when  required.  It  also  use  a  more  optimized  repartition  over 
the  cycle  of  the  updates,  as  given  by  Algorithm  1. 
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Figure  21:  Number  of  clusters  in  terms  of  time 

5.4  Notes  on  clustering  and  scheduling 

The  clustering  and  scheduling  algorithms  rely  on  some  parameters  that  must  be  set  properly 
in  order  that  our  adaptive  tracking  system  performs  well.  The  main  parameters  that  are 
concerned  with  the  performance  of  the  tracking  system  are: 

1.  the  Euclidian  distance  threshold  L  for  creating  clusters, 

2.  the  cluster  update  cycle  length  hc,  and 

3.  the  schedule  period  length  ^  of  the  time-slice  and  round-robin  scheduling  strategies. 

Note  that,  in  this  report,  hc  and  if:  are  equal.  For  the  results  presented  in  this  chapter, 
L,  hc  and  were  set  arbitrarily  to  provide  good  performance  with  the  three  scenarios. 
However,  in  practice  an  adaptive  tracking  system  would  require  L,  hc  and  il>  to  be  adjusted 
continually  and  based  on  the  situation  at  hand  (i.e.,  the  state  of  the  observed  targets). 
From  the  above  results,  we  already  know  that  ^  should  increase  as  the  number  of  targets 
increases,  such  that  all  targets  be  included  at  least  once  in  each  update  schedule.  Also,  L 
should  depend  on  the  targets’  distribution  and  on  the  sensor  characteristics.  For  example, 
if  for  a  particular  application  it  is  given  that  the  observed  targets  tend  to  be  very  close  in 
space  and  that  the  sensor  is  very  accurate,  then  L  should  be  set  low. 
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Figure  22:  Average  track-to-track  Mahalanobis  distance  ( AM ) 


Further  work  should  concentrate  on  replacing  the  above  presumptions  with  some  more 
theoretical  relations  involving  the  state  of  the  observed  space,  the  sensor(s)  characteristics 
and  the  required  performance. 
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Figure  23:  Total  number  of  Discrimination  Conflicts  ...,15]) 


Figure  24:  Average  determinant  of  error  covariance  matrix  ( ADET ) 
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6  Conclusion 


The  target  tracking  application  and  the  discrimination  power  problem  were  used  to  illus¬ 
trate  control  and  adaptation  concepts  in  data  fusion  applications.  A  two-loop  adaptation 
structure  was  defined  to  tackle  the  specific  adaptive  tracking  problem.  The  outer  loop  dy¬ 
namically  defines  and  manages  the  clusters  based  on  the  concept  of  discrimination  power, 
while  the  inner  loop  exploits  the  clusters  for  the  adaptation  purposes.  The  presented  re¬ 
sults  showed  the  superiority  of  the  time-slice  adaptation-based  strategy  over  static  policies 
such  as  the  round-robin  scheduling  strategy  that  was  tested  in  this  work.  Furthermore,  the 
time-slice  method  has  shown  to  be  superior  to  the  minimum  intra-cluster  discrimination 
power  method  since  the  latter  is  myopic  (i.e.,  does  not  plan  more  than  one  update  at  a 
time). 

Overall,  the  results  presented  in  this  report  showed  that  the  discrimination  power  can  be 
improved  by  adapting  the  target  tracking  operations.  In  the  light  of  these  results,  further 
works  involving  discrimination  power  should  concentrate  on.  The  following  list  gives  a 
tentative  ranking  of  problems  to  be  investigated.  This  ranking  is  based  on  the  challenges 
they  pose  and  the  preferences  of  the  authors. 

1.  Identify  and  study  specific  military  Command  and  Control  applications  on  which 
discrimination  power  might  have  an  impact.  For  example,  an  application  to  target 
engagement  in  dense  environments,  such  as  littoral,  is  being  developed.  The  objective 
is  to  measure  the  impact  of  an  improved  discrimination  power  on  the  frequency  and 
number  of  collateral  damages. 

2.  Extend  the  concept  of  discrimination  power  beyond  the  tracking  to  include,  for  in¬ 
stance,  the  identification  and  the  intent  and  capability  assessment. 

3.  Integrate  and  experiment  the  discrimination  power  control  strategies  on  some  typical 
surveillance  systems,  using  realistic  military  parameters  and  specifications. 

4.  Identify  and  study  the  inherent  factors  that  might  impact  on  the  discrimination  power. 

5.  Refine  the  clustering  and  scheduling  strategies  such  that  their  variable  parameters  be 
adjusted  depending  on  the  state  of  the  observed  space  and  on  the  state  and  charac¬ 
teristics  of  the  surveillance  system. 

6.  Integrate  the  discrimination  power  control  with  other  surveillance  and  sensing  objec¬ 
tives,  which  should  require  finding  new  sensor  scheduling  methods. 

7.  Refine  the  sensor  models  and  consider  various  types  of  sensors  working  together. 


DRDC  Valcartier  TR  2008-016 


35 


This  page  intentionally  left  blank. 


36 


DRDC  Valcartier  TR  2008-016 


List  of  symbols 


a  Degree  of  belief  (Bayesian  coverage  intervals)  or  significance 
level  (x2  hypothesis  testing) 
r)  Number  of  updates 
r  Process  gain  matrix 
A  Distance  threshold 
4>  Update  cycle 
v  Process  noise 
£  Constant 

(  Track  update  control  function 
ADET  Average  error  covariance  matrix  determinant 
AM  Average  Mahalanobis  distance 
C  Cluster  label 
d  Statistical  distance 
V  Dissimilarity  matrix 
DC  Discrimination  Conflict 
/  Frequency 

/(.)  State  transition  function 
F  Jacobian  of  f 
h  Track  update  period 
h(.)  Measurement  function 
hc  Outer-loop  period 
hs  Inner-loop  period 
H  Jacobian  of  h 
k  Discretized  time  instant 
K  Number  of  clusters 
L  Distance  threshold 
n  Loop  level  or  track  label 
nd  Number  of  dimensions  of  the  state  vector 


N  Number  of  tracked  targets  in  a  scenario 
P  Covariance  matrix  of  the  state  estimate 
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q  Process  noise  power  spectral  density 
Q  Track  quality 
R  Region  label 

R  Measurement  error  covariance  matrix 
s  Track  similarity  measure 
S  Schedule  vector 
t  Time  label 
w  Measurement  noise 
x  State  vector 
X  Abscissa 
x  State  estimate  vector 
Y  Ordinate 
z  Measurement  vector 
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List  of  acronyms 


c2 

CI 

EKF 

MICDP 

TP 

VOI 


Command  and  Control 
Coverage  Interval 
Extended  Kalman  Filter 

Minimum  Intra-Cluster  Discrimination  Power 
Tactical  Picture 
Volume  Of  Interest 
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Glossary 


The  followings  define  a  few  expressions  that  are  used  in  this  work.  Some  of  them  where 
inspired  from  [3,  31]. 

•  Accuracy 

—  Estimate  accuracy:  The  accuracy  associated  with  an  estimation  represents  its 
uncertainty  with  respect  to  the  true  value.  Under  the  normality  assumption,  it  is 
given  by  the  standard  deviation  of  the  estimate  posterior  probability  distribution 
function.  Equivalently,  accuracy  can  also  be  expressed  as  the  standard  deviation 
of  the  estimate’s  normally  distributed  random  error  with  mean  0.  Note  that 
accuracy  is  often  used  to  calculate  the  coverage  interval  of  estimates. 

—  Sensor  accuracy:  In  case  of  a  sensor,  the  accuracy  is  defined  according  to  the 
measurement  error  probability  distribution,  which  is  assumed  Gaussian.  Prac¬ 
tically,  the  value  of  accuracy  specified  in  most  sensing  systems  represents  the 
standard  deviation  of  the  distribution. 

•  Sensor  resolution:  The  resolution  of  a  sensor  is  its  ability  to  distinguish  between 
objects  that  are  very  close  in  space.  For  example,  a  typical  two-dimensional  (2D)  radar 
has  its  resolution  defined  for  both  range  and  bearing.  Accordingly,  range  resolution 
is  the  ability  of  a  radar  system  to  distinguish  between  two  or  more  objects  on  the 
same  bearing  but  at  different  ranges.  Angular  resolution  is  the  minimum  angular 
separation  at  which  two  equal  objects  can  be  separated  when  at  the  same  range.7 

•  Similarity:  Similarity  expresses  a  resemblance  between  two  or  more  estimates  with 
respect  to  the  estimated  features.  It  is  represented  by  statistical  measures  that  quan¬ 
tify  the  correlation  between  two  or  more  estimates.  For  example,  the  resemblance  can 
be  quantified  by  evaluating  the  amount  of  the  overlapping  portion  of  the  probability 
distributions  of  the  features. 

•  Dissimilarity:  Dissimilarity  is  the  opposite  of  similarity.  It  expresses  a  difference 
between  two  or  more  estimates  with  respect  to  the  estimated  features.  For  example, 
the  difference  can  be  quantified  by  evaluating  the  amount  of  the  non-overlapping 
portion  of  the  probability  distributions  of  the  features. 

•  Discrimination:  The  process  whereby  two  or  more  observed  objects  are  distin¬ 
guished. 

•  Discrimination  Power:  The  capability  (or  power)  to  distinguish  two  or  more  ob¬ 
served  objects.  In  target  tracking,  it  is  a  measure  of  the  dissimilarity  between  multiple 
target  tracks  distributed  in  the  space. 

7In  [38],  resolution  is  also  defined  as  the  finest  change  in  input  value. 
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